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ABSTRACT 

We present high-resolution HST images of all 35 AGNs with optical reverberation-mapping results, which 
we have modeled to create a nucleus-free image of each AGN host galaxy. From the nucleus-free images, we 
determine the host-galaxy contribution to ground-based spectroscopic luminosity measurements at A5100A. 
After correcting the luminosities of the AGNs for the contribution from starlight, we re-examine the H(3 R^i^g^-L 
relationship. Our best fit for the relationship gives a powerlaw slope of 0.52 with a range of 0.45-0.59 allowed 
by the uncertainties. This is consistent with our previous findings, and thus still consistent with the naive 
assumption that all AGNs are simply luminosity-scaled versions of each other. We discuss various consistency 
checks relating to the galaxy modeling and starlight contributions, as well as possible systematic errors in the 
current set of reverberation measurements from which we determine the form of the Rblr-L relationship. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: photometry — galaxies: Seyfert 



1. INTRODUCTION 

One of the key developments in extragalactic astronomy 
over the past decade has been the discovery that supermas- 
sive black holes are present in most, if not all, galaxies hav- 
ing a stellar bulge. Remarkably, the mass of the black hole 
is tightly correlated with the stellar velocity dispersion of the 
host g alaxy bulge (Ferrarese & Merritt 2000; Gebhard t et al.l 
I2OOOI) . pointing to a close link between the growth and evo- 
lution of galaxy stellar populations and the growth of nuclear 
black holes. Determining the masses of black holes in active 
galaxies is a crucial step toward understanding this connec- 
tion, and provides fundamental insight into the physics of ac- 
cretion and emission processes in the black hole environment. 

Unfortunately, most active galactic nuclei (AGNs) are too 
distant for black hole masses to be measured using spa- 
tially resolved stellar or gas dynamics. The technique 
that has been most successful for the measurement of the 
black hole mass (Mbh) in AGNs is reverb eration mapping 
(iBlandford & McKedfigsl lPetersonlfT993h . With this tech- 
nique, the AGN continuum (typically measured at 5 100 A) 
and broad emission lines (most notably, H/3) are monitored 
over an extended time period. Since the emission-line re- 
gions are photoionized by the central source, changes in the 
AGN continuum strength are followed by changes in the 
emission-line fluxes, with a time lag that depends on the 
light-travel time across the broad-line region (BLR). This 
time lag can be measured by cross-correlation of the con- 
tinuum and emission-line light curves, and gives the ra- 
dius of the BLR. Combining the BLR radius with the broad 
emission-line velocity width then gives the virial mass en- 
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close d within the BLR, wh i ch is domina ted by the black hole 
(e.g., iPeterson et all 1 19981: iKaspi et al.i r2000). The validity 
of reverberation masses has been upheld by the detection of 
virial behavior in the broad line r egion in a subs^ of ob- 
jects (e^gJPeterson & Wa ndelll999il2000HOnken & PetersonI 
l2002t lKollatschny..2003.) . as well as the consistency of rever- 
beration masses wit h other dynamical rnass methods, such 
as stellar dynamics ("DaviesetaT '200^, Onken et"an |2007|) 
and gas dynamics (Hicks & Malkan 2008|). Due to the long- 
term nature of reverberation-mapping projects, these mea- 
surements have only been carried out for a relatively small 
sample of AGNs in the past: about 36 Seyferts Is and low- 
luminosity quasars. 

The BLR radius-luminosity correlation {Rblr oc L") de- 
rived from this reverberation sample is the basis for all sec- 
ondary techniques used to est imate black hole masses in 
distan t AGNs (e.g.. lLaodll998t IWa ndel. Pet erson & MalkanI 
119991; iMcLure & JarvisI 120021: IVestergaard & PetersonI |2006|) 
and is an essential tool used to search for cosmologic al evo- 
lution of t he Mbh-c* relationship (e.g., " Peng et all 120061: 
IWoo et al.l I2OO8) . The power of the Rblr-L relationship 
comes from the simplicity of using it to quickly estimate Mbh 
for large samples of objects, even at high redshift, using two 
si mple measurem ents from a single spectrum of each object. 

iPeterson et al.l ("2004) compiled and consistently reanalyzed 
the database of available reverberation-mapping data for 35 
AGNs, with the goal of improving the measurements of the 
size of the BLR and th ereby impro ving their mass mea- 
surements. Subsequently, iKaspi et all ([2005) reanalyzed the 
Rblr-L relationship and found a power-law slope of a = 
0.665 ± 0.069 using the optical continuum and broad H/3 line. 
However, many of the AGNs in the sample reside in host 
galaxies that are comparable in luminosity to the AGN itself. 
Even worse, with the typically large apertures employed in 
reverberation-mapping campaigns (i.e. 5."0 x 7. "5), the host- 
galaxy starlight contribution to the spectroscopic luminosity 
of any given source is substantial. Failure to account for the 
enhancement of the luminosity by starlight results in an ar- 
tificially steep slope for the R^lr-L relationship, as a larger 
percentage of the luminosity in faint AGNs is contributed by 
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the host galaxy. 

A preliminary study by iBentz et aP (l2006ah presented 
two-dimensional fits to high-resolution HST images of 14 
reverberation-mapped AGNs, from which the host-galaxy 
contribution was determined. The luminosities of the 14 
sources were corrected and the R^ur-L relationship was re- 
examined, resulting in a measured power-law slope of a « 
0.5, consistent with the naive prediction that all AGNs are 
simply luminosity-scaled versions of each other In this work, 
we present high-resolution HST images of the rest of the H/3 
reverberation-mapped AGNs, bringing the total to 35. We 
improve upon our previous two-dimensional fits and reana- 
lyze the /^BLR-i- relationship after correcting the luminosity 
of every AGN for the contribution from starlight. We show 
that these new results are consistent with those from our pre- 
liminary study, and discuss the new measurements in light 
of known systematic errors that may affect the slope of the 
^BLR-i relationship. Throughout this work, we will assume a 
standard flat ACDM cosmology with ilfi = 0.04, rioivi = 0.26, 
1)a = 0.70, and Ho = 70km s"' Mpc"'. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Hubble Space Telescope 

Between 2003 August 22 and 2007 January 17, we ob- 
served 30 AGNs fro m the reverberation-mapped sample of 
iPeterson et al.l (l2004h with the HST Advanced Camera for 
Surveys (ACS). Following the failure of ACS on 2007 Jan- 
uary 27, the remaining 5 AGNs in our sample were observed 
with Wide Field Planetary Camera 2 (WFPC2). The targets 
are listed in Table 1 and details of the observations are listed 
in Table 2. 

For the ACS observations, each object was imaged with the 
High Resolution Channel (HRC) through the F550M filter 
(Ac = 5580 A and A A = 547 A), thereby probing the contin- 
uum while avoiding strong emission lines. The observations 
consisted of at least three exposures for each object, with ex- 
posure times of 120 s, 300 s, and 600 s. This method of gradu- 
ating the exposure times was employed to avoid saturation of 
the nucleus but still obtain a reasonable signal-to-noise ratio 
(S/N) for the wings of the point-spread function (PSF) and 
the host galaxy. Each individual exposure was split into two 
equal sub-exposures to facilitate the rejection of cosmic rays. 
The WFPC2 observations were centered on the PC chip and 
were taken through the F547M filter, the closest analog to the 
ACS F550M filter. The exposure times were again graduated, 
however, we used steps of 5s, 20s, 60s, 160s and 300s due 
to the smaller time-to-saturation afforded by the pixels in the 
PC chip. 

The data quality frames provided by the HST pipeline were 
consulted to identify the individual saturated pixels associated 
with the nucleus in each exposure frame. These saturated pix- 
els were clipped from the image and replaced by the same 
pixels from a non-saturated exposure after scaling them by 
the relative exposure times. All of the frames for each object 
were then summed to give one frame with an effective expo- 
sure time as listed in Table 2. 

Cosmic rays were identified in the summed images with 
the Laplacian cosmic ray identification package L. A. Cosmic 
dvan Dokkumll2001i) . Pixels in the PSF area of each image 
that were identified by L. A. Cosmic were excluded from the 
list of affected pixels prior to cleaning with XVista.^ Each 
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remaining affected pixel was replaced with the median value 
for the eight pixels immediately surrounding it. 

Finally, the summed, cleaned ACS images were corrected 
for the distortions of the camera using the PyRAF routine py- 
drizzle in the STSDAS^ package for IRAF The final stacked, 
cleaned images for all 35 AGNs are shown in Figure 1, over- 
laid with the spectroscopic aperture geometries from their 
ground-based monitoring campaigns. 

2.2. MDM Observatory 

Images of the reverberation-mapped sample of galaxies 
were also taken with the 1.3-m McGraw-Hill Telescope at 
MDM Observatory. Templeton, a 1024 x 1024 pixel CCD, 
was employed for the observations, giving a field-of-view 
(FOV)of8f53 X 8. '53 and a pixel scale of 0.50" pixel"'. Each 
galaxy was imaged through Harris B, V, and R filters. We fo- 
cus here on the observations of the seven NGC objects that 
were visible from the location of MDM Observatory^. A log 
of the MDM observations for those objects is presented in 
Table 3. The data were reduced and combined with IRAF'' 
following standard procedures. The final V-band images are 
shown in Figure 2. 

3. GALAXY DECOMPOSITIONS 

The images of each of the objects were modeled with typ- 
ical galaxy parameters in order to determine and accurately 
subtract the contribution from the central point source. The 
models were constructed usin g the two-dimens ional image 
decomposition program Galfit (iPeng et al.ll2002h . which fits 
analytic functions for the components of the galaxy, plus an 
additional point source for the nucleus, convolved with a user- 
supplied model point-spread function (PSF). 

For each of the objects in this study, the final cleaned, 
stacked image was fit with a central PSF and a constant sky 
contribution, as well as host- galaxy compo nents that were 
modeled using variations of the lSersid(ll968l) profile, 

I](r) = I],exp-'^["'/''''""-" (1) 

where rg is the effective radius of the component. Eg is the 
surface brightness at r^ n is the power-law index, and k, 
is coupled to n such that half of the total flux is within r^. 
Two special cases of the Sersic function are the exponen- 
tial pro file (n = I), often used in modeling galactic disks, 
and the Ide Vaucouleursi ([1948) profile (n = 4), historically 
used for modeling galactic bulges. We modeled disk com- 
ponents using the exponential profi le. Howev e r, we im- 
prove upon the results presented by iBentz et al] (l2006 a^ in 
that we employed the more general Sersic function for mod- 
eling bulges and we allow for additional parameters to de- 
scribe other surface brightness components (such as a bar 
or inner bulge component). These modifications were par- 
tially prompted by the large body of observations that find 
disk galaxies (of which our sample is mostly comprised) are 

maintained in the public domain by former Lick graduate students as a service 
to the community. It is currently maintained by Jon Holtzman at New Mexico 
State University, and is available at http://ganymede.nmsu.edu/holtz/xvista 
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more accurately describe d with bulge profiles that have n < 
4 (e.g. I^r mendv & Br uzuall 1 1 9781: IShaw & Gilmord 1 1 989t 
lAndredaids & Sandersi.l994l:lPeng et al.ll2002h . 

There is a paucity of archival stellar images with the HRC 
through the F550M filter, and so simulated PSFs were cre- 
ated using the Tiny Tim package (Krist 1993) which models 
the optics of HST plus the specifics of the camera and filter 
system. We tested our fits using a white dwarf image from the 
HST ai-chive as the PSF model (GO 10752, PI Lallo). Unfor- 
tunately, the image did not have the extremely high dynamic 
range necessary for a good fit when compared to these bright 
AGNs. 

Many of the nearest galaxies in this sample fill the FOV of 
the ACS HRC, and so the sky contribution could not simply 
be measured from the edges of the images. Rather, the sky 
level was iteratively determined, starting with an estimate of 
the sky brightness from the ACS Instrument Handbook as the 
initial input. The sky value was held fixed while additional 
parameters were fit to the galaxy. If the estimated sky value 
was too low, the Sersic index for the bulge would run up to the 
maximum allowed value, punching a "hole" in the nucleus 
of the image. If the estimated sky value was too high, the 
Sersic index could run down to zero, causing Galfit to crash. 
A sky value intermediate to these two situations was chosen 
such that the residuals and x~ values were minimized. Once a 
preliminary fit was achieved, the sky value was checked and 
adjusted if necessary, after which the galaxy parameters were 
re-fit. This is an improvement over our previous work, where 
the sky contribution was assumed to be negligible compared 
to the bright host galaxies of the 14 AGNs in that sample, 
especially as our expanded sample includes several AGNs that 
are significantly brighter than their host galaxies. 

All of the images required at least one galaxy component 
in addition to the sky and central PSF. Most required two, 
and a few required three or more to fit an additional compo- 
nent, such as a bar. Several of the bright AGNs required a 
small (< 2 pixel effective radius) component to help account 
for PSF mismatch between the AGN and the model PSF (for 
a full discus sion of PSF variations and mismatches in HST 
imaging, see lKrisjl2003l:lKlm et alll2008h . Extraneous objects 
in the field were also fit, such as intervening stars or galaxies, 
to ensure that the proper distribution of light was attributed 
to the AGN host galaxy. Compared to the simple galaxy de- 
compositions for the 14 objects in our previous work, the sur- 
face brightness residuals and nominal x~ fitting values have 
been significantly reduced, showing that the fits presented 
here more accurately model the underlying host galaxy sur- 
face brightness distributions. For example, the average mag- 
nitude of the residuals for 3C 120 was reduced by a factor 
of ~ 4 and the value was reduced by a factor of ^5.'" 
In all cases, the fits have been encouraged to attribute more 
flux to the sky background and PSF components, resulting in 
conservative flux values for the host-galaxy components. The 
quoted brightness for each of the host galaxy components may 
be somewhat underestimated as a result. 

For the seven NGC objects with MDM images, the V-band 
MDM images were each fit with an exponential disk and a 
Sersic bulge. The effective radius of the disk component was 

The actual value of i** not particularly meaningful when re- 
lated to galaxy fitting. The relative values for different fits to 
the same image are more relevant, with smaller values denoting "bet- 
ter" fits. Details of the calculation and interpretation of x" values 
can be found in the FA Q section o f Chien Peng's Galfit home page, 
|http://users. ociw.edu/peng/ work/galfit/ galfit.htmH 



then translated to the pixel scale of the HRC camera and held 
fixed during the fits to the HST image. As the seeing in the 
MDM images was typically on the order of ~ 2", the PSF and 
the bulge of the galaxy were blurred together Thus, the bulge 
fits from the ground-based images were held to be unreliable 
and were instead determined from the //ST images. 

Table 4 presents the details of the fits to the HST im- 
ages. The input image, Galfit model, residuals, and one- 
dimensional surface brightness cut for each of the 35 host 
galaxies are presented in Figure 3. Global parameters includ- 
ing the total galaxy luminosity and the ratio of the bulge lu- 
minosity to the total luminosity (fi/T) were determined from 
the fits and are listed in Table 5. Also listed in Table 5 are the 
morphological classifications of the galaxies, several of which 
were listed in the NASA/IPAC Extragalactic Database (NED). 
For those objects without morphological classifications (ob- 
jects marked with a flag in Table 5), the parameters fit to the 
galaxy images were us ed to determine t he ap propriate clas- 
sification based on the Ide VaucouleursI (11959 ) classification 
scheme. The index k for the subtype of the elliptical galaxies 
in the sample was calculated as k= \Q(\-b/d) and rounded 
to the nearest whole number. Spiral galaxies were classified 
based on their B/T compared to the mean of the distributions 
of BIT as a function of morphological type in Figure 6 of 
lKen3([T985h . 

As discussed by iPeng et all (l2002l) . degeneracy between 
galaxy components and between parameters within compo- 
nents is typically an issue when fitting analytic models to 
galaxy images. This is certainly the case with the sample of 
objects presented here. The bright AGNs in the galaxy centers 
are degenerate with the concentration of the bulge (the Ser- 
sic index n), especially when the bulge component is rather 
compact. As discussed above, many of these images have the 
added complication of having a somewhat uncertain sky con- 
tribution which also affects n as well as the scale length of the 
disk. Interpreting the specific details of the galaxy fits pre- 
sented here is therefore difficult. The relatively low n values 
for the galaxy bulges in this sample certainly seem to agree 
with the works referenced above that find n < 4 for most spi- 
ral galaxies. However, the actual n values may be somewhat 
underestimated as a result of our conservative fits that attribute 
more flux to the AGN and the sky background. A quick glance 
through the galaxy fits in Table 4 might lead one to specu- 
late on the prevalence of pseudobulge s versus classical bulges 
in this sample. Per the discussion bv lKormendv & KennicuttI 
(2004), however, the lack of dynamical info for the host galax- 
ies in this sample as well as the uncertainty in interpreting 
the galaxy fitting parameters complicates any conclusions that 
might be drawn about the origins of the bulges in these galax- 
ies. 

Several of the fits listed in Table 4 include components 
noted as being an "inner bulge" or "bar". The classifications 
for these extra components are simply morphological, with a 
"bar" typically having a more elongated shape than an "in- 
ner bulge". Bars themselves tend to have an inner bulbous 
component, so the "inner bulges" listed here may actually be 
related to bars themselves (Peng et al. 2002). We do not have 
kinematic information to determine which of these additional 
components may be physically distinct from the bulge or disk 
of the host galaxy. The fact that additional components are 
required to achieve a better fit is likely in many cases to sim- 
ply be the result of attempting to fit analytic functions to re- 
solved galaxies with strong substructure. In particular, areas 
with dust obscuration and spiral arm structure are difficult to 
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properly fit with analytic functions. As our main goal is to 
determine the PSF contribution as accurately as possible and 
to create "nucleus-free" images of these galaxies, we do not 
discuss here the origin or meaning of these additional compo- 
nents in the galaxy fits. 

4. FLUX MEASUREMENTS 

Once the fit to each galaxy was finalized, the sky and 
PSF components were subtracted, leaving a nucleus-free im- 
age of each host galaxy. Each image was overlaid with the 
typical aperture used in the ground-based monitoring pro- 
gram(s), at the typical orientation and centered on the posi- 
tion of the AGN (see Table 6). The counts within the aper- 
ture were summed and converted to f\ flux density units 
(ergs s"' cm-2 A"') using the keyword PHOTFLAM and 
the effective exposure time for each object. 

Color corrections to the observed galax y flux densities were 
calculated using a model bulge spectrum (iKinney et al.lll996h 
to account for the difference between the effective wavelength 
of the HST filter and restframe 5100 A for each object. The 
model bulge spectrum was redshifted to the distance of each 
AGN host galaxy and reddened by the appropriate Galactic 
extinction. The redshifted, reddened models were convolved 
with the HST filter response using the synphot package in the 
IRAF/STSDAS library and subsequently scaled to the appro- 
priate flux level as measured from the HST images. Finally, 
the observed flux at rest- frame 5 100 A was measured from the 
model. The flux density measured directly from the nucleus- 
free HST images and the color-corrected flux and luminosity 
at 5 100 A, are listed in Table 7. We note that the color correc- 
tion method outlined above is different from that employed 
previously, and again results in more conservative host galaxy 
flux values. 

The final step requires subtraction of the galaxy flux from 
the mean flux measured during a monitoring campaign, 
and leaves only the flux contribution from the AGN itself. 
The corresponding AGN luminosities were calculated and 
corrected for Galactic absorption using the ISchlegel et al.l 
([T998i) An values lis ted in NED and the extinction curve of 
ICardelh et all (1 1989b . adjusted to Av/E{B-V) = 3.1. The 
starlight-corrected AGN fluxes and luminosities are listed in 
Table 8 with their corresponding H/3 lags. The numbers in 
bold font are the weighted averages of multiple measurements 
for a particular object. 

Ongoing work with this sample of objects has resulted 
in some up dates to the reverbe ration datab a se pre sented by 
iPeterson et al. (2004) and fit bv lKaspi etall (|2005|) . In addi- 
tion, a few data sets require cautionary or explanatory notes. 
Specifically, we note the following: 

• NGC3516: The light curves presented by 
IWanders et al.l (Il993h were not measured from 
data with an absolute flux calibration. Therefore, while 
we have included this object, it should be regarded 
with caution as the spectroscopic luminosity is not 
well-determined. 

• NGC4151: The ape rture used i n the m onitoring cam- 
paign described by iMaoz et al.l (1199 Ih is 20" x 28", 
which is larger than t he field of vi e w of t he HRC cam- 
era. In addition, the iKaspi et aH (1 19961) data set has 
a rather unconstrained time lag that is consistent with 
zero once the monotonic increase in the contin uum and 
line flux is removed (see lMetzroth et aLll2006l for a full 



discussion). Both of these data sets for NG C4151 have 
been s uperseded by the results reported by iBentz et aP 
(I2006bt) . which we include here in their stead. 

• PG1211H-143: IPeterson et aH (|2004 note that afl of 
the lag measurements for this object are rather uncon- 
strained. For consistency with past analyses of the 
R^hR-L relationship, we have included the H/3 lag for 
PG121 1 here, but it should be regarded with caution. 

• NGC4593: T Tie ca mpaign originally described by 
iDietrich et al.l (I1994J) and reanalyzed by ' Onken et al.l 
(l2003h gives an H/3 lag that is consistent with zero. The 
results for NGC 4593 have bee n superseded by those re- 
ported by Penney et al.l (l2006l) . which we include here. 

• IC 4329A: According tol Peterson etaLl(l2004 . the light 
curves for this object are very noisy and of poor quality, 
resulting in an H/3 lag measurement that is highly sus- 
pect as well as being consistent with zero. In addition, 
IC 4329A is an edge-on galaxy with a prominent dust 
disk along the line of sight to the AGN, resulting in a 
substantial amount of internal reddening. The presence 
of the strong dust lane across the galaxy in the HST im- 
ages also presented problems for the host galaxy fitting. 
Because this object has both an unreliable lag measure- 
ment and an unreliable luminosity measurement, we 
have excluded it from the analysis of the R^ug^-L re- 
lationship. 

• Mrk279: The m onitoring campaign described by 
iMaoz et al.1 (Il990h used an aperture of 20" x 28", 
which exceeds the HRC field of view. We in- 
clude here only the results for Mrk279 reported by 
Hantos-Lleo et al. (2001). 

• NGC 5548: The apertur e of 20" x 28" used in the cam- 
paign described by Netz er et al.l(ll990h is larger than the 
field of view of the HRC camera. We i nclude here the 
14 ind ivid ual measurements r eported by IPeterson et al.l 
(I2002h and lBentzetan(l2007h . 

• PG2130H-099: For some time now, the H/3 lag mea- 
sured for this object has been suspected of being ar- 
tificially long. Th is suspicion has been confirmed by 
iGrier et aH (l2008l) . who have recently reanalyzed the 
data set originally presented by lKaspi et al.l (1 19961) and 
find that the lag measured from these data is erroneous 
due to undersampling of the light curves and long-term 
secular changes in the equivalen t width of the broad 
H/3 emission line. In addition, Grier et all present a 
new monitoring campaign for this object that was un- 
dertaken at MDM Observatory in the fall of 2006, and 
report a new H/3 lag measurement. We include only the 
iGrier et aTl result in our analysis here. 

5. CONSISTENCY CHECKS 

The difficulty of fitting surface brightness profiles with mul- 
tiple analytic functions is widely appreciated. While the 
GALFIT program we have used works remarkably well and 
represents the state of the art, there are still some ambiguities 
that cannot always be easily resolved: at a sufficiently large 
distance, for example, it is simply impossible, even with the 
high angular resolution of the ACS HRC, to separate a point 
source AGN and a luminous bulge with a high Sersic index. 
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In addition, several parameters that are used in the surface 
brightness fits can be very degenerate, leading to difficulties in 
disentangling the results from many different parameter com- 
binations when no other information is available. All this may 
result in systematic uncertainties that are larger than the for- 
mal errors (of order 10%) that are listed in Table 7. There are, 
however, a number of consistency checks that provide robust 
upper limits on the host galaxy contributions and can help to 
remove some of this degen eracy. Indeed, our decision to re- 
visit the lBentz et al.l (l2006al) fits was motivated at least in part 
by such considerations. 

The simplest consistency check is that the AGN flux must 
always remain non-negative: the host galaxy flux at any wave- 
length must be less than the observed brightness of the com- 
bined AGN and host galaxy when the AGN is in its faintest 
observed state. We have compared the measured host galaxy 
fluxes to the individual monitoring light curves for each of the 
objects in this sample to ensure that this consistency check is 
always satisfied. A second consistency check is provided by 
spectral decomposition, i.e., by fitting multiple components 
of a spectrum rather than an image. Our experience, however, 
is that spectral decomposition is certainly no less ambiguous 
than image decomposition, but when carefully done with very 
good data, shows reasonable consistency (e.g., the modeling 
of the optical specti'um of NGC5548 bv iDennev et al.ll2008l 
yields a host galaxy flux in agreement with the value we find 
here). A third check, which applies at least in the case of the 
higher luminosity AGNs, is that the total luminosity of the 
host galaxy, based on the model fit, must not exceed the lumi- 
nosity of known bright normal galaxies. Our galaxy fits seem 
to be consistent with this constraint as well. The brightest host 
galaxy in the sample is PG 1226H-023, also known as 3C273, 
which has My ~ -23.8, and the second brightest with My « 
-23.2 is PG 1700H-518, the mos t distant AGN in the s ample 
at z = 0.292. Compared to the iTrentham et"an ( 12005! ) field 
galaxy luminosity function, the derived host galaxy bright- 
nesses for the AGNs in this study seem to be fairly typical. 

A more subtle consistency check on the host-galaxy con- 
tribution to the ground-based monitoring luminosity comes 
from photoionization considerations, namely that the ratio of 
emission-line to AGN continuum flux (i.e., the emission-line 
equivalent width) ought not to vary wildly in time. For an 
optically thick medium, the equivalent width of an emission 
line is expected to stay relatively constant. Some variations 
are possible, either due to changes in the continuum shape or 
if there is an optically thin BLR component. Significant vari- 
ations in the equivalent width can follow a large-amplitude 
change in the incident continuum flux, however, the varia- 
tions should stabilize and return to a relatively constant value 
shortly thereafter. In addition, secular trends in the equiva- 
lent width can occur over dynamical timescales, which are 
much longer than reverberation timescales. If too much flux 
is attributed to the host galaxy of an object, the AGN contin- 
uum will be underestimated and the fractional variations of 
the H/9 emission-line equivalent width will become consis- 
tently large, in contrast to the variations that can be expected 
from the discussion above. An important point to note is that 
the relevant equivalent width measure is the ratio of emission- 
line flux at some time, f, relative to the continuum flux at 
some earlier time, f - r, where r is the emi ssion-line time de- 
lay from reverberation measure ments (e.g.. lPogge & PetersonI 
[1991 iGilbert & Peter son 2003^. 

To illustrate this check, we consider the specific case of 
3C 120, which we selected because the starlight contribution 



changed the most from our previous work, decreasing by a 
factor of almost 3. In Figure 4 we show the equivalent width 
of the Hp emission line (in the observ ed frame) as a functio n 
of time, based on the light curve from [Peterson et akl ( Il998l) . 
Each measurement of the emission-line flux, Fupit), was di- 
vided by an interpolated value of the continuum flux den- 
sity, F\{t - r), with an observed-frame emission-line lag of 
T = 39.4 days. We note in passing that the first few emission- 
line measurements are disregarded since there is no corre- 
sponding continuum information. The open circles in Fig- 
ure 4 refer to the original data, without any correction for 
starlight. The triangles are the new equivalent width mea- 
surements for the same emission line fluxes, but in this case 
the continuum measu rements have bee n adjusted by the host- 
galaxy value from B entz et al.l ('2006a). The large equivalent 
widths and consistently rapid variations strongly suggest that 
this value of the host-galaxy correction is too large, as the 
fractional variations of the remaining AGN continuum must 
be enormous. The filled circles show the equivalent widths 
based on the revised value of the host-galaxy contribution 
given here in column (3) of Table 7. In this case, the val- 
ues of the equivalent widths are much more reasonable (i.e., 
typical of quasars, where host-galaxy contamination is neg- 
ligible) and the variations are much more moderate, indicat- 
ing that the host-galaxy measurement provided here is more 
accurate than the measurement based on th e more simplistic 
host-galaxy models utilized by lBentz et al.l 

For all 35 objects in the sample, the galaxy fits in §3 were 
carried out independently from the consistency checks dis- 
cussed above. The more conservative approach we have ap- 
plied to the galaxy fits here results in good agreement with 
the expectations from all of these consistency checks for the 
35 objects in the sample. Of the 14 objects included in our 
preliminary study, three objects have starlight contributions 
through their monitoring apertures that changed by more than 
35% as measured directly in counts from the nucleus-free im- 
ages, one of them being 3C 120, as discussed above. In all 
cases, the decreased residuals and values from the new fits, 
combined with the results of the consistency checks, lead us 
to believe that the galaxy fits in this work more accurately 
describe the underlying host-galaxy surface bri ghtness dis- 
tributions than the simplistic fits determined bv iBentz et al.l 
(■2006a ). 

6. THE RADIUS-LUMINOSITY RELATIONSHIP 

We have calculated several fits to the Rqi^r-L relation- 
ship for the full sample of starlight-corrected AGNs. For 
this analysis, the time lags have been restricted to the H/3 
lag only. This is different from the method employed by 
iKaspi et akl (l2005h . where the /^blr-^' fits were also calcu- 
lated for the Balmer-line average, which sometimes included 
Ha and/or H7. There is one exception to this, however, in that 
PG0804H-761 does not have a reliable H/3 measurement, but 
it does have reliable measurements for Ha and H7. We use 
the lag measurement for Ha here, as it is the more reliable of 
the two measurement s available. 

In the same way as iKaspi et akl we have made the distinc- 
tion of treating each separate H/3 lag measurement of an ob- 
ject individually, as well as taking the mean of multiple mea- 
surements weighted by the average of the positive and nega- 
tive errors. We tested the differences between weighting mea- 
surements by the average of their errors, by taking only the 
positive e rrors, and by takin g the errors toward the fit in the 
manner of lKaspi et al.l ( l2005h . We find the differences in these 
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weighting methods to be at the 2% level, and therefore negli- 
gible. 

While the two above methods for sampling the 
reverberation-mapping database are relatively straight- 
forward to carry out, various issues arise when determing 
the RsLR-L relationship from these datasets. In the first 
case, where every measurement is given equal weight, those 
objects that have many measurements (such as NGC5548) 
will have more weight in the determination of the slope 
than objects with single measurements. In the second case, 
where multiple measurements are averaged together, each 
object has the same weight but information is being lost in 
the average because we do not, in fact, expect the luminosity 
and lag to be the same in multiple campaigns during different 
years. It is unclear as to the correct way to combine multiple 
measurements in this case. 

A somewhat more laborious method of sampling the 
reverberation-mapping database is to randomly select one pair 
of radius and luminosity measurements for every object, from 
which selection the relationship is fit, and to build up a large 
number (A^ = 1000) of individual realizations through Monte 
Carlo techniques. This method gives equal weight to every 
object and circumvents the problem of how to combine multi- 
ple measurements for any particular object. In addition to the 
two simpler methods outlined above, we employ this method 
in the determination of the Rblr-L relationship, but we con- 
sider this method to be superior to the others. 

Three different fitting routines were used to calculate the 
relationship between the size of the BLR and the optical lu- 
minosity: 

1. FITEXY (IPress e t aLlll992l) . which estimates the parame- 

ters of a straight-line fit through the data including er- 
rors in both coordinates. FITEXY numerically solves 
for the minimum orthogonal using an interative 
root-findi ng algorithm . We include intrinsic scatter 
similar to Kaspi et alJ (l2005l) after the prescription of 
[Tremaine et al. (200^. Namely, the fractional scatter 
listed in Table 9 is the fraction of the measurement 
value of Rblr (not the error value) that is added in 
quadrature to the error value so as to obtain a reduced 
of 1.0. 

2. BCES ( lAkritas & Bershadvl[T996 l). which attempts to ac- 

count for the effects of errors on both coordinates in the 
fit using bivariate correlated errors, including a com- 
ponent of intrinsic scatter. We adopt the bootstrap of 
the bisector value following iKaspi et aTl (l2005h . with 
A^= 1000 iterations." 

3. GaussFit jMcArthur et al.l 19941) . which implements gener- 

alized least-squares u sing robust Househo lder Orthog- 
onal Transformations ( Jeffervs 1980, 1981) to solve the 
non-linear equations of condition for the problem of er- 
rors in both coordinates. No attempt is made to account 
for intrinsic scatter. 

Table 9 lists the fit parameters determined for a R^lr-L rela- 
tionship of the following form: 

log(/?BLR) = ^+alog(ALA(5100A)) (2) 

' ' We note that the errors on L listed in the table are formal measurement 
errors that, in most cases, are far smaller than the actual observed variations 
even in one campaign. Therefore, we also tested BCES correlations that as- 
sume an uncertainty of 0. 1 -0. 1 5 dex on L. The results are similar to the ones 
obtained by our realization method and hence not listed separately. 



where a is the slope of the power-law relationship between 
Rblr and AL(5100A), and K is the zero point. The calcu- 
lated power-law slopes to the R^lr-L relationship range from 
0.499 ± 0.042 to 0.554 ± 0.050, depending on the particular 
alorithm used for fitting the relationship and how the objects 
with multiple measurements are treated. 

We prefer the Monte Carlo random sampling method out- 
lined above as the proper way to treat objects with multi- 
ple individual measurements. And for its manner of dealing 
with intrinsic scatter within the data set, we prefer the BCES 
bootstrap method. The combination of these methods gives 
K = -21.3!|^ and a = 0.519!!J:;!^^, our best estimate for the 
form of the Rblr-L relationship at this time. 

Figure 5 shows the Rblr-L relationship after correcting 
the full sample of reverberation-mapped AGNs for the con- 
tribution from host-galaxy starlight. The top panel of Fig- 
ure 5 shows each individual data point from the monitoring 
campaigns included here, and the bottom panel shows a sin- 
gle data point for each individual AGN determined by the 
weighted average of multiple measurements. The solid lines 
show the best fit to the Rblr-L relationship described above. 

7. DISCUSSION 

Throughout this work, we have improved upon our original 
methods by using more accurate and more conservative pro- 
files to model the host galaxy components, as well as a more 
conservative color correction method, all resulting in conser- 
vative measurements of the host galaxy starlight as measured 
for every object contributing to the Rblr-L relationship. Even 
so, the best fit to the relation ship has no t changed signifi- 
cantly from that presented by iBentz et al.1 : 0.5 19;!^ com- 
pared to the previous value of 0.5 18 ±0.039. It would appear 
that given reasonably high-resolution, unsaturated images of 
AGN host galaxies and reasonable fits to the host galaxy sur- 
face brightness profiles, the AGN luminosity can be corrected 
for the contamination from starlight fairly accurately. 

AGN BLRs can be modeled assuming the central radia- 
tion field is the only source of heating and ionization of the 
gas. The simplest, most naive assumption is that all BLRs 
are made of identical clouds, with the same density, column 
density, composition and ionization parameter. In addition, 
the radiation from all AGNs can be assumed to have the same 
spectral energy distribution (SED). This results in the predic- 
tion that Rblr oc L"'^. Any changes in the BLR gas distribu- 
tion or the AGN SED, especially those associated with source 
luminosity, will result in a different slope. In fact, AGNs 
have been obs erved to have different SEDs as a function of lu- 
minosity (e.g..'Mushotzky & Wandell ll989l:IZheng & MalkanI 
IT993). It is therefore interesting that our present work based 
on the results from several reverberation-mapping campaigns 
produces a slope for the Rblr-L relationship that is consistent 
with the naively expected slope of 0.5, and that, to a first ap- 
proximation, brighter AGNs are simply "scaled up" versions 
of fainter AGNs. 

Removing the host-galaxy starlight component reduces the 
scatter'^ in the Rblr-L relationship, fr om 39-44% as fou nd 
by the best-fit FITEXY results from iKaspi et al.1 (l2005h to 
34-40% as found here. It also flattens the slope of the re- 
lationship considerably. This has the overall effect of bias- 
ing samples that used previous versions of the Rblr-L re- 
lationship to estimate black hole masses. The host-galaxy 

Throughout this manuscript, "scatter" refers to the I— a deviation from 
the best fit Sblr-^- relationship. 
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starlight is often removed through spectroscopic decomposi- 
tion of single-epoch spectra before the size of the BLR ra- 
dius is estimated from the continuum luminosity. However, 
it is crucial to use a Rslr-L relationship in which the objects 
providing the calibration are also corrected for host-galaxy 
starlight, which we h ave provided here. Compared with the 
best-fit relationship of lKaspi et alj(l2005h . we find that the cal- 
ibration of the Rbi^r-L relationship presented here results in 
black hole masses that are 30% smaller at L = 10'**', 50% 
larger at L = 10'*'*, and a factor of ^ 3 larger at L = lO'*^. 

At the low-luminosity end (L < 10"*^), there may still be 
some uncertainty as to the behavior of the Rblr-L relation- 
ship. In the current sample of objects with reveberation- 
mapping results, there are somewhat fewer objects at lower- 
luminosity and they have larger uncertainties than the other, 
higher-luminosity objects in the sample. It should be kept 
in mind that the lower-luminosity objects were, in general, 
the first targets of ground-based monitoring campaigns due 
to their relatively low redshifts and high apparent bright- 
nesses. The larger uncertainties in their measurements is par- 
tially due to the less-rigorous control over observational fac- 
tors in those early monitoring campaigns (such as observing 
cadence, spectral resolution, detector efficiency, etc.) simply 
because there was a lack of experience in this field at that 
time. The problems of small sample size and relatively larger 
uncertainties for lower luminosity objects will soon be miti- 
gated by two independent reverberation mapping campaigns 
which have been recently carried out at MDM Observatory 
and at Lick Observatory and targeted the low-luminosity end 
of the relationship. Preliminary results from the MDM cam- 
paign promise to replace several m easurements, as was done 
in the case ofNGC 4593 dPennev e t al. 2006) and NGC4151 
dBentz et al.ll2006bl) in the 2005 MDM campaign. And pre- 
liminary results from the Lick campaign promise to add sev- 
eral new objects to the low-luminosity end of the Rslr-L re- 
lationship (e.g., Bentz et al. 2008). 

Internal reddening is known to be a problem in some of 
the very nearest, and lower-luminosity, objects in the current 
reverberation sample. For example, NGC 3227 is known to 
have substantial internal reddening compared to most of the 
other objects in the current sample and sho uld therefore re- 
quire o ne of the largest reddening corrections. ICrenshaw et al.l 
(1200 Ih determined a reddening curve for NGC 3227 and find 
thatat5100A,AA/£'(B-V) = 3.6 and ^(B-y) = 0.18. This 
implies an extinction at 5100 A of 0.65 mag, which, if cor- 
rected for, would increase the luminosity here by a factor of 
1 .8, or 0.26 dex, moving the location of NGC 3227 in the bot- 
tom panel of Figure 5 from slightly left of the best fit R^lr-L 
relationship to right on top of it. As we do not have simi- 
lar corrections for the other objects in this sample, we do not 
apply the reddening correction for NGC 3227 in our deter- 
mination of the Rblr-L relationship. However, based on the 
magnitude of the correction determined for NGC 3227, there 
does not seem to be any reason to expect that correcting all 
the sources for internal reddening will have much effect on 
the slope of the R^lr-L relationship. 

An obvious issue that may also be addressed with the 
galaxy fits that we have presented in this work is the relation- 
ship between black hole mass and host galaxy bulge luminos- 
ity (or mass; 'K ormendy & Richstondfl995t iMagorrian et al.l 
[1998.) . We discuss this rela tionship for the AG Ns in our sam- 
ple in a related manuscript (iBentz et al.ll2009l) . 



We have presented high-resolution HST images of the 35 
AGNs with optical reverberation-mapping results. The host 
galaxy of each object was fit with typical galaxy compo- 
nents, and a nucleus-free image of each AGN host galaxy 
was created. From these nucleus-free images, we measured 
the starlight contribution to the ground-based spectroscopic 
luminosity measured at rest- frame 5100 A. We then removed 
the starlight contamination from the AGN luminosities and 
re-examined the Rbi^r-L relationship. We find a best fit slope 
of a = 0.519^Q Qgg, consistent with the results from our pre- 
liminary study, and still suggesting that all AGNs are simply 
luminosity-scaled versions of each other. We discuss several 
consistency checks that support our galaxy modeling results. 
Various systematics, such as the smaller number and larger 
uncertainties of measurements at lower luminosities as well 
as internal reddening, are discussed in the context of their ef- 
fect on the Rblr-L relationship. 
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Fig. 1. — The final stacked images for the full sample of 35 reverberation-mapped AGNs with /fSr imaging. The ground-based spectroscopic monitoring 
apertures are overlaid. Each image is 25" x 25", with north up and east to the left. 
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Fig. 2. — V-band images of the seven NGC objects in tlie reverberation-mapped sample that are visible from MDM Observatory. Each image is 5' x 5', with 
north up and east to the left. 
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Fig. 3. — From left to right, each row shows the following: the HST image of an AGN and its host galaxy; the best-fit to the galaxy+AGN surface brightness 
profiles from Galfit; the residuals of the fit; isophotal analysis of the image and models, with the data points measured from the sky-subtracted HST image, the 
solid line from the total host-galaxy model image (which is convolved with the PSF for comparison with the observed image), and the dashed line from the PSF 
model for the AGN component. 
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Fig. 3. — Continued. 
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Fig. 3. — Continued. 
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Fig. 4. — Equivalent width consistency clieck for 3C 120. The equivalent width of the H/3 emission line is shown as a function of time for three cases: with 
no correction of the continuum flux density for host-galaxy contamination (open circles); with corrections for the host galaxy using the value from 'Bentz et al.l 
(2(306a) (triangles); and with corrections for the host galaxy using the revised value presented here in Table 7 (filled circles). The larger host-galaxy correction 
from B entz et al.l leads to equivalents widths that are unacceptably large and variable. The new, more conservative value for the starlight contribution corrects 
this problem so that the equivalent width values and variations are in keeping with expectations from photoionization physics. 
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AL^ (5100 A) 

Fig. 5. — The H/3 Rblr-L relationship after correcting the AGN luminosities for the contribution from host-galaxy starlight. The top panel shows each separate 
measurement as a single data point, and the bottom panel shows the weighted mean of multiple measurements for any individual object. The solid lines are the 
best fit to the relationship (listed in bold face in Table 9), which has a slope of a = 0.519;^ q^^. 
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\J\J 1 J .\J 


-1-11 IfiOl 

T 1 ,J 1 U J 


14900 
yj. i'tzuw 


671 7 


307 




pn 00^2+2^1 


on 54 52 1 


+25 25 3R 

1 -jCj 


1 5500 


740.0 


0.205 




H^iiTnll Q 
1 all till y 


01 9^ 4S 780 


-58 4S 90 50 


04709 


208 6 


1 16 




Mrk SQO 
iviiiv ^y\j 




-00 Af, 00 00 


09(il0 


115 3 


161 


Nnn 8(1 1 

ly V 1 V ou J 


3C 120 


04. 1 1 OOSS 


-1-05 91 15 (i90 


03301 


145 


1 283 


Mrif 1 50(1 

iVlllv 


Arlc 190 


OS 1(1 1 1 491 


-00 08 50 18 


03271 


141 8 


554 


Mrk 1005 
ivii iv i.yjy^ 


Mrk 79 


07 49 ^9 707 


-1-4.0 48 14. 75 


09910 


96 7 


305 




pn 0S04.+7/il 


08 10 'SS flOO 


-1-7^1 09 49 00 

T / U \JZ "tZ . 17 VJ 


1 0000 


460 5 


150 






08 47 42 4 


-1-14.45 04 


0(i400 


287 4 


159 




Mrk 1 1 

iVlllS. 1 IVj 


OQ 95 19 S70 


-1-59 17 1 59 

TJZ 1 / lKJ.JZ. 


01590 

JZ7 


155 2 


056 




pn oo^'^+did 


OQ S(i 59 4 


+41 15 22 


23410 


1 172 1 


054 




Nnr "^997 


10 9^ ^0 57Q0 


-1-10 51 54 1 80 

T17 »J1 J'T.lOl.7 


0018(1 


23.6 


098 






11 0fi47 400 


-1-79 14 C\f\ 88 
T /Z JM- VJO.OO 


00884 


JO. 1 


183 






11 ni 79 


17 44 1 8 Q 


n 00071 




514 
U. J i'4 




Nnr 4.0 s 1 


1 9 0^ OQ fil4 


-1-4.4 11 59 80 

T't't J 1 J Z . Ow 


00914 

U.WUZ J'T 


15.2 


056 




Nnr 4. 1 s 1 


19 1 019 S7Q 


+39 24 20 63 

TJ^ Zt- ZVJ.UJ 


00119 


14.3 


1 19 




pn 191 i+M'^ 

1 W IZ-llTlT-J 


12 14 17 7 


-1- 1 4. 01 19 (i 

TlT-UJ IZ.U 


08000 
yj.yjoyjyyj 


367.6 


150 




pn 1 99(1-1-09 


1 9 90 Ofi ^iOQ7 


4-09 01 08 508 


15834 


757 5 


089 


3C 273 


pn 1 99Q-I-904. 


1 9 19 0'^ ^i05 

iZ JZ U-J.UVjJ 


-1-90 00 90 91 

TZU U7 Z?.Z 1 


0(ilO 1 


282 8 


1 17 


Mrk 771 Mr Tnn 1 549 

iVlllS. Ill LX. HJIl U'tZ 


i N V 1 V J ^ J 


1 9 10 10 495 

iz jy j!!7.'tzj 


-05 90 10 14. 

vjJ zvj jy.jH- 


00000 


38 8 


106 


Mrk 1 110 

iVll JV 1 J J u 


pn 1 "^07-1-08 s 


1 1 00 47 

1 J \jy / .w 


-i-OS 10 480 


1 5500 


739 2 


145 




IC 4329A 


1 1 40 10 9(1 

1 -J 'T7 1 7.ZU 


-10 1 8 14 

J\J lO J'T.'J 


1 (i05 


69.6 


255 




Mrk 279 


1 1 51 01 447 


-i-fiO 1 8 90 57 

10Z7.J/ 


01045 
yj.yjjyj^^ 


133 5 


068 




PG 141 1+442 


141348.3 


+440014 


0.08960 


409.7 


0.036 




NGC5548 


1417 59.534 


+25 08 12.44 


0.01718 


74.5 


0.088 




PG 1426+015 


14 29 06.588 


+01 17 06.48 


0.08647 


394.4 


0.137 




Mi-k817 


14 3622.068 


+5847 39.38 


0.03146 


137.9 


0.029 


PG 1434+590 


PG 1613+658 


16 13 57.179 


+65 43 09.58 


0.12900 


605.6 


0.114 


Mrk 876 


PG 1617+175 


16 2011.288 


+17 2427.70 


0.11244 


521.9 


0.180 


Mrk 877 


PG 1700+518 


17 0124.800 


+514920.00 


0.29200 


1509.6 


0.151 




3C 390.3 


184208.9899 


+794617.127 


0.05610 


250.5 


0.308 




Mrk 509 


204409.738 


-1043 24.54 


0.03440 


151.2 


0.248 




PG 2 130+099 


2132 27.813 


+10 08 19.46 


0.06298 


282.6 


0.192 


IIZw 136 & Mrk 1513 


NGC 7469 


23 03 15.623 


+08 5226.39 


0.01632 


70.8 


0.297 


Mrk 1514 



^ Dis tances were calculated from the redshifts of the objects, except for NGC 3227 - where we use the distance measured by surface brightness fluctuations to NGC 3226 1 Blakeslee et al] 
with which NGC 3227 is interacting — and NGC 405 1 — where we use the average TuUy-Fisher distance reported bv lRusselli i i20Q3t) .^ Values are from iSchlegel et alj |199S) . 
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TABLE 2 
HST Observation Log 



Object 


Observational 


Date Observed 


Beginning UTC 


Total Exposure 


Dataset 




Setup 


(yyyy-mm-dd) 


(hh:mm) 


Time (s) 




IVliJv J J J 


.1.1 1 1 W - 1 JJV^iVl 


900f^-0S-94 


08 "26 


2040 


TQ1VTTT010 




WFPr9 F^47M 

VV 1 JT Jt^ / iVl 


9007-Ofi-Ofi 


21-39 


1445 


TTQMTT590 


IT V J J 1 1 


AnS T-TRn FSSOM 


900fi-0S-9S 


09:59 


2040 


TQMT T030 


H^llTnll 

1 £1110.11 y 


AnS HRn F^SOM 


90O'^-OS-99 


00:44 


1020 


T8SP040 


Mrlr son 




9001 19 IS 


02"27 


1020 




3C 120 


XCK HRn F^SOM 


900'^-19-OS 


05 "48 


1 020 


T8*\P0(S0 


AVn 190 

rt.lvll 1 Z.\.> 


An*^ HRn FSSOM 


900^1-10-10 


18"36 


2040 


TQMT T540 


Mrk79 




900fS 1 1 OS 


18"33 


2040 


TQ1VTTT050 


pn nKnd+7fi i 


AnS T-TRn FSSOM 


900fS-0Q-90 


23:09 


2040 


TQMT T060 




An<5 HRn FSSOM 


9004. OS 10 


20*11 


1090 


TR<2P100 


Mrlf 1 1 




9004-0^-98 


17:34 


1020 


TRSP1 10 




ACS HRr FSIflM 


2006-10-2*5 


18:51 


2040 


T9MTT070 




APS HRP FSSflM 


9004-0'^-90 


04:28 


1020 


J OkJ^-.l jU 


i> VJV_. J J lU 


APS HRP FSSnM 


900^-1 9-1 Q 


01:03 


2220 


TQDOOIO 


Mnr ■^TS'^ 


rtv^o ,ni\.^^ ,r J J Uivi 


9001 11 \'\ 




1090 


TRSPl '^0 






9004-09-16 


01:49 


1020 


TSSP1 60 


Nnr 41^1 


APS HRP FSSOM 




14-25 


1020 


TSSP1 70 


pn 191 1 + 1 4'^ 


AP'^ HRP FSSOM 

,^\. .J.I 1 IW - 1 jjvyivi 




02"30 


2040 


TQMT TOSO 


L vJ 1 Z.Z.UTVZ. J 


APS HRP FS'SOM 

.7.1 IIW -I jjvyivi 


9007-01-17 


12-40 


2040 


TQMT TOQO 


pn 1 990-1-904. 


APS HRP FSSOM 


900/^ 1 1 90 


02'41 


2040 


TQMT T 1 00 




APS HRP F^SOM 


900^1-01 -10 


21-05 


2220 


TQnn090 


pn 1 "^OT-i-oss 


WFPP9 F^47M 

VV 1 IT / iVl 


9007-01-91 


14:36 


1445 


TT0MTT1 10 




APS HRP F^SOM 


900fi-09-99 


00:03 


2220 


TQnno40 


Mrk 279 


ACS,HRC,F550M 


2003-12-07 


03:54 


1020 


J8SC240 


PG 141 1+442 


ACS,HRC,F550M 


2006-11-10 


23:52 


2040 


J9MU120 


NGC5548 


ACS,HRC,F550M 


2004^04-07 


01:53 


1020 


J8SC270 


PG 1426+015 


WFPC2,F547M 


2007-03-20 


16:18 


1445 


U9MU130 


Mrk817 


ACS,HRC,F550M 


2003-12-08 


18:08 


1020 


J8SC290 


PG 1613+658 


ACS,HRC,F550M 


2006-11-12 


04:44 


2040 


J9MU140 


PG1617 175 


WFPC2,F547M 


2007-03-19 


17:59 


1445 


U9MU150 


PG 1700+518 


ACS,HRC,F550M 


2006-11-16 


07:51 


2040 


J9MU160 


3C 390.3 


ACS,HRC,F550M 


2004-03-31 


06:56 


1020 


J8SC340 


Mrk 509 


WFPC2,F547M 


2007-04-01 


22:50 


1445 


U9MU170 


PG2130+099 


ACS,HRC,F550M 


2003-10-21 


06:47 


1020 


J8SC360 


NGC7469 


ACS,HRC,F550M 


2006-07-09 


22:00 


2220 


J9DQ030 



TABLE 3 
MDM Observation Log 













Exposure Times 






Object 


Date Observed 


Airmass 


B 


V 


R 


Seeing 




(yyyy-mm-dd) 


(sec 


z) 


(s) 


(s) 


(s) 


(") 


NGC3227 


2003-02-08 


1.023- 


1.135 


1800 


1680 


795 


1.77 


NGC3516 


2003-04-25 


1.317- 


1.392 


1900 


845 


1500 


1.85 


NGC4051 


2003-02-08 


1.048- 


1.111 


1250 


795 


690 


1.71 


NGC4151 


2003-02-09 


1.008- 


1.050 


1470 


1200 


1220 


2.07 


NGC4593 


2003-02-09 


1.298- 


1.576 


1650 


1860 


1380 


2.58 


NGC5548 


2003-04-25 


1.018- 


1.308 


4000 


2500 


1080 


1.63 


NGC 7469 


2003-09-26 


1.087- 


1,106 


1440 


1 260 


1 560 


1 .69 
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TABLE 4 
Details of Galaxy Fits 



Object 


Sky 
(counts) 


a 


i?e 




n 


b/a 


Note 


Mrk335 


10.0 


14.9 












PSF 






16.3 


0.51 


6 


.61 


0, 


,84 


Bulge 






15.8 


1.49 


1 


.00 


0. 


,95 


Disk 


PG 0026+129 


2.8 


16.1 












PSF 






17.0 


0.16 


2, 


.58 


0. 


,52 


Add'l PSF 






17.3 


7.21 


1 


.72 


0. 


,86 


Bulge 


PG 0052+251 


5.7 


15.6 












PSF 






18.0 


0.17 


3, 


.12 


0. 


,34 


Bulge 






17.2 


8.24 


1 


.00 


0. 


,69 


Disk 






18.8 


3.77 


5, 


.47 


0, 


78 


Field galaxy 


Fairall9 


7.0 


15.1 












PSF 






15.3 


0.49 


5 


.61 


0. 


,94 


Bulge 






15.3 


2.92 


1 


.00 


0. 


,44 


Disk 


Mrk590 


6.8 


17.9 












PSF 






16.3 


0.44 


1 


.22 


0. 


,62 


Inner bulge 






15.8 


0.75 





.59 


0. 


,96 


Bulge 






14.2 


3.14 


1 


.00 


0. 


,91 


Disk 


3C120 


8.0 


14.9 












PSF 






18.3 


0.04 


3 


.67 


0. 


10 


Add'l PSF 






17.6 


0.66 


1 


.10 


0. 


,89 


Bulge 






16.1 


3.11 


1 


.00 


0, 


,65 


Disk 


Ark 120 


12.3 


14.7 












PSF 






15.0 


0.05 


3, 


.62 


0. 


,84 


Bulge 






14.9 


1.85 


1 


.00 


0. 


,87 


Disk 


Mrk79 


13.7 


15.7 












PSF 






16.1 


0.85 


2, 


.79 


0. 


,67 


Bulge 






16.6 


1.98 


1 


.00 


0, 


66 


Disk 






14.7 


8.89 





.55 


0. 


,24 


Bar 


PG 0804+761 


10.7 


15.0 












PSF 






14.7 


0.08 


1 


.22 


0. 


,73 


Add'l PSF 






16.8 


3.33 


1 


.00 


0, 


,74 


Bulge 


PG 0844+349 


6.0 


14.6 












PSF 






16.9 


0.04 


2 


.28 


0. 


12 


Bulge 






16.7 


2.86 


1 


.00 


0. 


,75 


Disk 


MrkllO 


3.7 


16.1 












PSF 






18.2 


0.25 


1 


.35 


0, 


,85 


Bulge 






16.5 


1.73 


1 


.00 


0. 


,93 


Disk 






16.6 












Star 


NGC3227 


14.2 


15.2 












PSF 






15.4 


0.01 


1 


.51 


0, 


,92 


Add'l PSF 






14.9 


0.06 


1 


.08 


0, 


,68 


Inner bulge 






12.8 


1.19 


2, 


.14 


0. 


,50 


Bulge 






13.8 


4.66 


1 


.00 


0. 


,47 


Disk 


NGC3516 


12.0 


15.2 












PSF 






13.4 


0.38 


1 


.24 


0, 


,77 


Inner bulge 






13.0 


1.74 





.96 


0, 


60 


Bulge 






14.4 


4.53 


1 


.00 


0. 


,52 


Disk 


PG 0953+414 


10.0 


14.9 












PSF 






17.7 


28.82 


1 


.39 


0. 


,55 


Bulge 


NGC3783 


12.0 


14.2 












PSF 






14.7 


0.49 


1 


.09 


0. 


,92 


Bulge 






15.0 


1.95 





.33 


0. 


,29 


Bar 






12.0 


6.02 


1 


.00 


0. 


,83 


Disk 


NGC4051 


8.0 


14.8 












PSF 






15.1 


0.03 


1 


.07 


0. 


,85 


Inner bulge 






15.1 


0.07 





.31 


0. 


,71 


Inner bulge 






12.8 


0.86 


1 


.80 


0, 


,49 


Bulge 






12.4 


4.24 


1 


.00 


0. 


,67 


Disk 


NGC4151 


0.2 


14.5 












PSF 






14.4 


0.07 


4, 


.29 


0. 


,54 


Inner bulge 






14.0 


0.14 





.71 


0. 


,96 


Inner bulge 






12.0 


0.73 





.81 


0. 


,95 


Bulge 






13.0 


3.77 


1 


.00 


0, 


69 


Disk 


PG 1211+143 


17.0 


22.5 












PSF 






14.9 


0.06 





.08 


0. 


,63 


Add'l PSF 






17.3 


2.59 


1 


.00 


0. 


,84 


Bulge 


PG 1226+023 


6.5 


15.4 












PSF 






13.2 


0.14 





.22 


0, 


,29 


Add'l PSF 






15.6 


11.71 


2, 


.50 


0. 


,75 


Bulge 


PG 1229+204 


15.6 


16.7 












PSF 






20.2 


0.23 





.82 


0. 


,69 


Bar? 






17.3 


0.99 


1 


.15 


0, 


,87 


Bulge 






15.8 


7.84 


1 


.00 


0. 


,62 


Disk 


NGC4593 


15.0 


15.3 












PSF 






15.1 


0.52 





.09 


0. 


71 


Inner bulge 
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UDject 


(counts) 


nistmag^ 


(kpc) 


n 


U In 

u[a 


Note 






12.3 


2.83 


1.94 


0.68 


Bulge 






13.5 


9.68 


1.00 


0.52 


Disk 


PG 1307+085 


1.3 


15.6 








PSF 






18.1 


0.19 


0.43 


0.05 


Add'l PSF 






17.6 


8.74 


1.25 


0.80 


Bulge 


IC4329A 


15.0 


14.6 








PSF 






13.8 


0.01 


0.31 


0.87 


Add'l PSF 






13.7 


0.77 


0.39 


0.96 


Inner bulge 






12.6 


3.36 


0.50 


0.43 


Bulge 






15.8 


8.74 


1.00 


0.41 


Disk 


Mrk279 


6.0 


15.0 








PSF 






16.2 


0.05 


5.94 


0.65 


Inner bulge 






16.2 


0.97 


1.88 


0.60 


Bulge 






15.2 


3.14 


1.00 


0.54 


Disk 


PG 1411+442 


9.3 


15.7 








PSF 






15.8 


0.05 


1.00 


0.48 


Add'l PSF 






16.8 


5.05 


1.71 


0.58 


Bulge 






18.8 


1.87 


2.23 


0.64 


Field galaxy 


NGC5548 


0.2 


16.7 








PSF 






14.7 


1.18 


4.36 


0.86 


Inner bulge 






13.8 


2.95 


1.39 


0.90 


Bulge 






15.5 


11.51 


1.00 


0.85 


Disk 


PG 1426+015 


1.8 


15.2 








PSF 






16.3 


4.74 


1.63 


0.77 


Bulge 






19.8 


0.91 


3.14 


0.56 


Field galaxy 






21.7 


0.15 


1.98 


0.50 


Field galaxy 


MrkSl? 


4.7 


15.1 








PSF 






17.7 


0.28 


2.44 


0.81 


Bulge 






14.4 


4.32 


1.00 


0.74 


Disk 


PG 1613+658 


11.0 


15.2 








PSF 






16.3 


6.75 


1.35 


0.80 


Bulge 






19.1 


1.79 


3.64 


0.60 


Field galaxy 


PG 1617+175 


1.3 


16.3 








PSF 






17.4 


0.15 


0.09 


0.07 


Add'l PSF 






17.3 


3.28 


5.35 


0.84 


Bulge 


PG 1700+518 


9.0 


19.2 








PSF 






15.3 


0.15 


0.02 


0.58 


Add'l PSF 






16.9 


0.28 


0.37 


0.65 


Add'l PSF 






17.7 


12.06 


5.61 


0.89 


Bulge 


3C 390.3 


1.3 


15.8 








PSF 






17.0 


0.77 


3.86 


0.74 


Bulge 






16.9 


2.54 


1.00 


0.86 


Disk 


Mrk509 


4.2 


14.5 








PSF 






15.4 


0.04 


0.03 


0.48 


Add'l PSF 






15.0 


1.85 


1.00 


0.79 


Bulge 


PG 2130+099 


4.0 


14.8 








PSF 






18.9 


0.38 


0.56 


0.37 


Bulge 






16.5 


4.66 


1.00 


0.55 


Disk 


NGC7469 


7.0 


15.2 








PSF 






15.0 


0.37 


1.37 


0.70 


Inner bulge 






13.4 


3.68 


1.31 


0.55 


Bulge 






14.6 


6.99 


1.00 


0.94 


Disk 






16.2 








Star 



The magnitude of an object is computed as: ms,„„j = -2.5 log (^^^) + zpt, vihexezpt = lAASl for the F550M filter with ACS HRC andzpf = 21.685 for the F547M filter with the 
PC chip on WFPC2. 
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TABLES 
Global Galaxy Parameters 



Object ALxaaiCSlOO)" B/T Morphological Flag"^ 

(10"'* ergs s"') i 2^ Classification 



Mrk335 


0.37 


0.37 




so/a 




PG 0026+129 


2.3 


1.00 




El 


* 


PG 0052+251 


4.1 


0.32 




Sb 




Fairall 9 


2.4 


0.52 




SBa 


* 


Mrk590 


1.4 


0.17 


0.28 


SA(s)a 




3C120 


0.73 


0.21 




SO 




Aknl20 


2.1 


0.49 




Sbpec 




Mrk79 


0.73 


0.19 


0.87 


SBb 




PG 0804+761 


1.5 


1.00 




E3 


* 


PG 0844+349 


1.2 


0.45 




Sa 


* 


MrkllO 


0.26 


0.17 




Sc 


* 


PG 0953+414 


3.6 


1.00 




E4 


* 


NGC3227 


0.23 


0.65 


0.74 


SAB(s) pec 




NGC3516 


0.69 


0.52 


0.86 


(R)SB(s) 




NGC3783 


1.5 


0.07 


0.12 


(R')SB(r)a 




NGC4051 


0.16 


0.36 


0.45 


SAB(rs)bc 




NGC4151 


0.20 


0.60 


0.75 


(R')SAB(rs)ab 




PG 1211+143 


0.59 


1.00 




E2 


* 


PG 1226+023 


11.3 


1.00 




E3 


* 


PG 1229+204 


1.8 


0.20 


0.21 


SBc 


* 


NGC4593 


0.91 


0.71 


0.76 


(R)SB(rs)b 




PG 1307+085 


1.8 


1.00 




E2 


* 


IC4329A 


2.4 


0.70 


0.96 


SAO 




Mrk279 


0.95 


0.21 


0.43 


SO 




PG 1411+442 


1.2 


1.00 




E4 


* 


NGC5548 


0.98 


0.60 


0.87 


(R')SA(s)0/a 




PG 1426+015 


1.8 


1.00 




E2 


* 


Mrk817 


1.2 


0.05 




SBc 




PG 1613+658 


4.2 


1.00 




E2 


* 


PG 1617+175 


1.3 


1.00 




E2 


* 


PG 1700+518 


6.7 


1.00 




El 


* 


3C 390.3 


0.87 


0.48 




Sa 


* 


Mrk509 


0.95 


1.00 




E2 


* 


PG 2130+099 


0.88 


0.10 




(R)Sa 




NGC7469 


1.4 


0.64 


0.79 


(R')SAB(rs)a 





Galaxy luminosities are determined directly from the model parameters that were fit to the HST images. They do not include features that were not modeled, such as the nuclear 
starburst ring in NGC 7469, and do not include the contributions from field galaxies or stars. Bulge luminosities here include the contribution from any bar and/or inner bulge 
component.'^ Morphological classifications are from NED where available. Those marked with a flag are determined from the galaxy fit parameters in this work, as described in the 
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TABLE 6 

Ground-Based Monitoring Apertures and Observed Fluxes 



Object Reference" PA Aperture /((l+z)5100) 

(°) ( " X ") (10-'5 ergs cm'^ A"' ) 



Mrk335 


1 


90.0 


5.0 


X 


7.6 


7.68 ±0.53 




1 


90.0 


5.0 


X 


7.6 


8.81 ±0.47 


PG 0026+ 129 


2 


42.0 


10.0 


X 


13.0 


2.69 ±0.40 


PG 0052+251 


2 


153.4 


10.0 


X 


13.0 


2.07 ±0.37 


Fairall 9 


3 


0.0 


4.0 


X 


9.0 


5.95 ±0.66 


Mrk590 


1 


90.0 


5.0 


X 


7.6 


7.89 ±0.62 




1 


90.0 


5.0 


X 


7.6 


5.33 ±0.56 




1 


90.0 


5.0 


X 


7.6 


6.37 ±0.45 




1 


90.0 


5.0 


X 


7.6 


8.43 ±1.30 


3C120 


1 


90.0 


5.0 


X 


7.6 


4.30 ±0.77 


Aknl20 


1 


90.0 


5.0 


X 


7.6 


10.37 ±0.46 




1 


90.0 


5.0 


X 


7.6 


7.82 ±0.83 


Mrk79 


1 


90.0 


5.0 


X 


7.6 


6.96 ±0.67 




1 


90.0 


5.0 


X 


7.6 


8.49 ±0.86 




1 


90.0 


5.0 


X 


7.6 


7.40 ±0.72 


PG 0804+761 


2 


315.6 


10.0 


X 


13.0 


5.48 ±1.00 


PG 0844+349 


2 


36.8 


10.0 


X 


13.0 


3.71 ±0.38 


MrkllO 


1 


90.0 


5.0 


X 


7.6 


3.45 ±0.36 




1 


90.0 


5.0 


X 


7.6 


3.96±0.51 




1 


90.0 


5.0 


X 


7.6 


2.64 ±0.86 


PG 0953+414 


2 


31.7 


10.0 


X 


13.0 


1.56±0.21 


NGC3227 


4 


25.0 


1.5 


X 


4.0 


23.46 ±3.70 




5 


90.0 


5.0 


X 


10.0 


12.70 ±0.68 


NGC3516 


6 


25.0 


1.5 


X 


2.0 


7.83 ±2.35 


NGC3783 


7 


0.0 


5.0 


X 


10.0 


11.38 ±0.95 


NGC4051 


8 


90.0 


5.0 


X 


7.5 


13.38±0.92 


NGC4151 


9 


90.0 


5.0 


X 


12.75 


23.8 ±3.0 


PG 1211+143 


2 


352.2 


10.0 


X 


13.0 


5.66 ±0.92 


PG 1226+023 


2 


171.2 


10.0 


X 


13.0 


21.30±2.60 


PG 1229+204 


2 


291.5 


10.0 


X 


13.0 


2.15 ±0.23 


NGC4593 


10 


90.0 


5.0 


X 


12.75 


15.9±0.7 


PG 1307+085 


2 


186.5 


10.0 


X 


13.0 


1.79±0.18 


IC4329A 


11 


90.0 


5.0 


X 


10.0 


5.79 ±0.73 


Mrk279 


12 


90.0 


5.0 


X 


7.5 


6.90 ±0.69 


PG 1411+442 


2 


347.0 


10.0 


X 


13.0 


3.71 ±0.32 


NGC 5548 


13 


90.0 


5.0 


X 


7.5 


9.92 ±1.26 




13 


90.0 


5.0 


X 


7.5 


7.25 ±1.00 




13 


90.0 


5.0 


X 


7.5 


9.40 ±0.93 




13 


90.0 


5.0 


X 


7.5 


6.72 ±1.17 




13 


90.0 


5.0 


X 


7.5 


9.06 ±0.86 




13 


90.0 


5.0 


X 


7.5 


9.76±1.10 




13 


90.0 


5.0 


X 


7.5 


12.09 ± 1.00 




13 


90.0 


5.0 


X 


7.5 


10.56± 1.64 




13 


90.0 


5.0 


X 


7.5 


8.12±0.91 




13 


90.0 


5.0 


X 


7.5 


13.47 ± 1.45 




13 


90.0 


5.0 


X 


7.5 


11.83±1.82 




13 


90.0 


5.0 


X 


7.5 


6.98 ±1.20 




13 


90.0 


5.0 


X 


7.5 


7.03 ±0.86 




14 


90.0 


5.0 


X 


12.75 


6.63 ±0.36 


PG 1426+015 


2 


341.4 


10.0 


X 


13.0 


4.62 ±0.71 


Mrk817 


1 


90.0 


5.0 


X 


7.6 


6.10±0.83 




1 


90.0 


5.0 


X 


7.6 


5.00 ±0.49 




1 


90.0 


5.0 


X 


7.6 


5.01 ±0.27 


PG 1613+658 


2 


164.2 


10.0 


X 


13.0 


3.49 ±0.43 


PG 1617+175 


2 


253.0 


10.0 


X 


13.0 


1.44 ±0.25 


PG 1700+518 


2 


183.5 


10.0 


X 


13.0 


2.20±0.15 


3C 390.3 


15 


90.0 


5.0 


X 


7.5 


1.73 ±0.28 


Mrk509 


1 


90.0 


5.0 


X 


7.6 


10.92 ± 1.99 


PG 2 130+099 


16 


0.0 


3.0 


X 


6.97 


4.63 ±0.23 


NGC 7469 


17 


90.0 


5.0 


X 


7.5 


13.57 ±0.61 



REFERENCES. — I JPeterson etay ^ j 19980 . 2. iKaspi et a i] j200CI) . 3. ISantos-Lleo et alj 11991. 4. ISalamanca _eta lJ 11994. 5. IWinge et alj I1995tl . 6. IWanders et al, 09931) . 7. 
IStirpeet alj 119941) . 8. Peterson et al. (2000), 9. Bentz et al. 1 2006b), 10. Denney et al. 1 2006). II. Winge et al. i mS). 12. ISantos-Lleo et alj 120011) . 13. iPeterson et alj 120021) and 
references therein, 14. Bentz et al. ( 2007), 15. Dietrich et al. ( 1998), 16. Grier et al. (2008), 17. . Collier et al.,a998i) . 

Note. — Here, and throughout, observed galaxy fluxes are tabulated at rest-frame 5100 A. 
^ References refer to reverberation-mapping campaigns in optical wavelengths. 
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TABLE? 

Host-Galaxy Fluxes and Luminosities 



Object 




fg^l(HST) 


/gal((l+Z)5100) 


AL;,.g,,(5100)'' 




(10"'^^ ergs s"' cm"2 A"') 


f-.i((l +7)5100) 


(10~^^ ergs s~^ cm""^ A~') 


(10'''' ergs s"') 


Mrk335 


1.88 


0.85 


1.60 ±0.15 


0.142 ±0.013 


PG 0026+129 


0.381 


1.00 


0.379 ±0.035 


1.46±0.13 


PG 0052+251 


0.713 


0.98 


0.699 ± 0.064 


3.08 ±0.28 


Fairall 9 


3.47 


0.88 


3.07 ±0.28 


0.927 ±0.085 


Mrk590 


4.81 


0.85 


4. 10 ±0.38 


0.384 ±0.035 


3C120 


0.783 


0.82 


0.641 ±0.059 


0.217 ±0.020 


Akn 120 


6.70 


0.85 


5.68 ±0.52 


1.079 ±0.099 


Mrk79 


1.74 


0.84 


1.46 ±0.13 


0.106 ±0.010 


PG 0804+761 


0.703 


0.97 


0.683 ±0.063 


1.076 ±0.099 


PG 0844+349 


1.24 


0.92 


1.14±0.11 


0.684 ±0.063 


MrkllO 


0.786 


0.88 


0.688 ±0.063 


0.109 ±0.010 


PG 0953+414 


0.208 


1.11 


0.231 ±0.021 


2.47 ±0.23 


NGC3227'' 


4.12 


0.82 


3.37 ±0.31 


0.0124 ±0.0011 




8.57 


0.82 


7.01 ±0.65 


0.0258 ± 0.0024 


NGC3516 


4.55 


0.82 


3.73 ±0.34 


0.0382 ±0.0035 


NGC3783 


6.06 


0.80 


4.86 ±0.45 


0.0776 ± 0.0071 


NGC4051 


10.1 


0.82 


8.28 ±0.76 


0.0123 ±0.0011 


NGC4151 


21.6 


0.82 


17.6±1.6 


0.0240 ± 0.0022 


PG 1211+143 


0.633 


0.95 


0.598 ±0.055 


0.592 ±0.054 


PG 1226+023 


1.37 


0.98 


1.34±0.12 


5.741 ±0.529 


PG 1229+204 


1.48 


0.92 


1.36±0.13 


0.766 ±0.071 


NGC4593 


10.7 


0.82 


8.85±0.82 


0.0889 ±0.0082 


PG 1307+085 


0.232 


1.00 


0.233 ±0.021 


0.986 ±0.091 


IC4329A 


4.43 


0.83 


3.67 ±0.34 


0.133 ±0.012 


Mrk279 


3.49 


0.87 


3.02 ±0.28 


0.355 ±0.033 


PG 1411+442 


0.826 


0.96 


0.791 ±0.073 


0.904 ± 0.083 


NGC 5548'' 


4.63 


0.84 


3.88 ±0.36 


0.143 ±0.013 




5.51 


0.84 


4.61 ±0.43 


0.169±0.016 


PG 1426+015 


1.19 


0.96 


1.14±0.11 


1.29±0.12 


Mrk817 


1.77 


0.87 


1.54±0.14 


0.188±0.017 


PG 1613+658 


1.52 


0.98 


1.50±0.14 


4.08 ±0.38 


PG 1617+175 


0.341 


0.99 


0.336 ±0.031 


0.701 ±0.065 


PG 1700+518 


0.246 


1.41 


0.347 ±0.032 


6.79 ±0.63 


3C 390.3 


0.945 


0.90 


0.853 ±0.079 


0.430 ±0.040 


Mrk509 


2.74 


0.92 


2.52 ±0.23 


0.435 ±0.040 


PG2130+099 


0.440 


0.92 


0.405 ±0.037 


0.240 ±0.022 


NGC 7469 


10.2 


0.83 


8.43 ±0.78 


0.327 ±0.030 



^ The galaxy luminosities presented here are measured through the ground-based monitoring aperture directly from the PSF-subtracted HST images. Any field galaxies or stars, or 
additional unmodeled galaxy structures, that are included in the original ground-based spectroscopic aperture contribute to this luminosity.'' The two different entries for NGC 3227 
and NGC 5548 correspond to the two different monitoring apertures that were employed during the spectroscopic monitoring programs for these objects. They are listed in the same 
order as in Table 6. 
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TABLE 8 

Rest-frame Time Lags and Starlight Corrected Luminosities 



Object 


H/3 Time Lag /agn((1 + 


s)5100A) ALa.agn(5100A) 




(days) (10^" ergs s~ 


'cm-^A-') (lO'^^ergss-') 



Mrk335 



PG 0026+129 
PG 0052+251 

Fairall 9 
Mrk590 



3C120 
Ark 120 



Mrk79 



PG 0804+761 
PG 0844+349" 
MrkllO 



PG 0953+414 
NGC3227 



NGC3516 
NGC3783 
NGC4051 
NGC4151 
PG 1211+143 
PG 1226+032 
PG 1229+204 
NGC4593 
PG 1307+085 
IC4329A'' 
Mrk279 
PG 141 1+442 
NGC5548 



PG 1426+015 
Mrk817 



PG 1613+658 
PG 1617+175 
PG 1700+518 
3C 390.3 



12.531 
15,7+3.4 

lll-O-i'l 
89.8!2« 

2o.rf 

29 2+*'^ 
28.8;3-« 
25.6+2" 



38.1+ 



^-15.3 
47^+8.3 

31. M 

39.7+1-9 

9.0!«| 
16.1+«-° 



16.0 
15.2 



146.9+11 



32.3![ 



24.3 
20.4: 
33.3+ 



9 

10.0 
8.2+5„l 

6.7;«-8 

10.2+3 
5.8!f 

6.6!J: 

93.8^:^5 

306.8^^ 

3.7+« 
105.6+3° 



37. i 



I. 51: 

16.7« 
124.3+^[ 
19.7![ 
18.6+1 
15.9+1 

II. 0+'' 



-2.0 

i3.o;} ' 

13.4«: 

21.7;| 

16.4:', 

17.5+f* 

26.5!^' 

24.8;f 

6.5:| 
14.3«' 

6.3+1 
18.0+" 



95.0 



1+2' 

-37, 

19.0:1 
15.3:| 
33.6!° 
21.8!^ 
40.1![f 
71.5+1 
251.8!*| 

23.6 



6.09 ±0.53 
7.21 ±0.47 
6.72 ±0.35 

2.31 ±0.40 
1.37 ±0.37 
2.88 ±0.66 
3.80 ±0.62 
1.23 ±0.56 
2.27 ±0.45 
4.3 ±1.3 
2.44 ±0.30 
3.66 ±0.77 

4.69 ±0.46 
2. 14 ±0.83 

4.09 ±0.40 
5.50 ±0.67 
7.03 ±0.86 

5.94 ±0.72 
6.03 ±0.43 

4.8±1.0 
2.57 ±0.38 
2.77 ±0.36 
3.28±0.51 

1.95 ±0.86 

2.83 ±0.28 
1.33±0.21 
20.1 ±3.7 

5.70 ±0.68 
6.17 ±0.67 

4.1 ±2.3 
6.52 ±0.95 
5.10±0.92 

6.2±3.0 
5.06 ±0.92 
20.0 ±2.6 
0.79 ±0.23 
7.05 ±0.70 
1.56±0.18 
2. 12 ±0.73 
3.88 ±0.69 
2.92 ±0.32 
6.0±1.3 
3.4±1.0 
5.52 ±0.93 
2.8±1.2 
5.19±0.86 
5.9±1.1 
8.2±1.0 
6.7 ±1.6 
4.25 ±0.91 
9.6 ±1.5 
8.0±1.8 
3.1 ±1.2 
3.16±0.86 
2.02 ±0.36 

3.84 ±0.23 
3.48 ±0.71 
4.56 ±0.83 

3.46 ±0.49 

3.47 ±0.27 
3.54 ±0.23 
1.99 ±0.43 

1.10 ±0.25 

1.85 ±0.15 
0.88 ±0.28 



0.541 ±0.047 
0.640 ±0.041 
0.603 ±0.031 

8.9 ±1.6 
6.0 ±1.6 

0.87 ±0.20 
0.355 ±0.058 
0.1 15 ±0.053 
0.212 ±0.043 

0.41 ±0.12 
0.287 ±0.032 

1.24 ±0.26 
0.889 ±0.088 

0.41 ±0.16 
0.847 ±0.081 
0.401 ±0.049 
0.513 ±0.063 
0.435 ±0.053 
0.447 ±0.031 
7.6± 1.6 

1.54±0.23 
0.439 ±0.058 
0.520 ± 0.080 

0.31 ±0.14 
0.461 ±0.045 

14.2 ±2.2 
0.074 ±0.014 

0.0209 ± 0.0025 
0.0304 ± 0.0031 
0.042 ±0.024 
0.104±0.015 
0.0076 ± 0.0014 
0.0084 ±0.0041 

5.00 ±0.91 
86±11 

0.45 ±0.13 
0.0708 ± 0.0070 

6.59 ±0.77 
0.077 ±0.026 
0.456 ±0.082 

3.34 ±0.36 
0.222 ±0.047 
0.124 ±0.037 
0.203 ±0.034 
0.105 ±0.043 
0.191 ±0.032 
0.216 ±0.040 
0.302 ±0.037 
0.246 ±0.060 
0.156 ±0.033 
0.352 ±0.054 
0.292 ±0.067 
0.1 14 ±0.044 
0.1 16 ±0.032 
0.074 ±0.013 
0.205 ±0.011 

3.94±0.81 

0.56±0.10 
0.423 ±0.060 
0.424 ± 0.032 
0.438 ±0.028 
5.4 ±1.2 

2.30 ±0.51 

36.3 ±2.9 
0.44 ±0.14 
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TABLES — Continued 



Object 


H/3 Time Lag 
(days) 


/agn((1+z)5100A) 
(10"'5 ergs s"' cmr^ A"') 


ALaagn (5 100 A) 
(lO'^'ergss-') 


Mrk509 
PG 2 130+099 
NGC 7469 


79.6+«' 
no q+4.7 

^■■^-0.8 


8.4±2.0 
4.22 ±0.23 
5.14±0.61 


1.45 ±0.34 
2.51 ±0.14 
0.200 ±0.023 



Note. — Numbers in boldface are the weighted averages of all the measurements for that particular object. Fluxes are as observed. Time lags and luminosities are listed in the 
rest-frame of the object with weighted averages calculated in log space. 

^ The H/3 lag measurement for this object was deemed unreliable by lPeterson et alj |2004j) . In its place, we give the most reliably measured lag, which is for Ha. Because of the 
extremely poor lag measurement (which is consistent with zero) and the poor flux calibration for this object, we do not include it in the fit to the Rblr-L relationship. 
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TABLE 9 
H/3RBLR-iFlTS 



Note N K a Scatter^ 



FITEXY 



All 59 -21.0±1.8 0.511 ±0.041 34.0 
Avg 34 -22.1 ±2.3 0.535 ±0.051 40.0 
MC 34 -22.3 ±2.2 0.540!«;j54 40.3!^i 



BCES 



All 59 -20.4±1.8 0.499±0.042 
Avg 34 -21.5 ±2.1 0.524 ±0.046 
MC'' 34 -21.3;|'^ 0-519!ljo^^ 



GaussFit 



All 59 -21.7±1.5 0.529±0.033 
Avg 34 -21.8±1.9 0.531 ±0.042 
MC 34 -22.9 ±2.2 0.554^0°^^ 



Note. — All: Each individual measurement is treated separately. Avg: Multiple measurements for a single source are combined into a weighted average. MC: Monte Carlo 
techniques are used to randomly sample the multiple measurements for a single source, producing one pair of Rbui and L measurements per object. The values and uncertainties 
presented for the fit using this method describe the median and 68% confidence intervals for the distributions of slopes and intercepts built up over multiple realizations. As described 
in the text, IC 4329A was not included in any of the fits listed here. 

^ The scatter listed here is the percentage of the measurement value of /?blr that is added in quadrature to the error value so as to obtain a reduced of 1 .0.'^ This fit, which properly 
treats multiple measurements of individual objects and accounts for intrinsic scatter in the data set, should be taken as our current best estimate for the form of the ^blr-^^ relationship. 



